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Abstract 

The validity of fluctuation dissipation relations in an aging system is studied 
in a colloidal glass during the transition from a fluid-like to a solid-like state. 
The evolution of the rheological and electrical properties is analyzed in the range 
IHz — 4:0Hz. It is found that at the beginning of the transition the fluctu- 
ation dissipation relation is strongly violated in electrical measurements. The 
amplitude and the persistence time of this violation are decreasing functions 
of frequency. At the lowest frequencies of the measuring range it persists for 
times which are about 5% of the time needed to form the colloidal glass. This 
phenomenology is quite close to the recent theoretical predictions done for the 
violation of the fluctuation dissipation relation in glassy systems. In contrast 
in the rheological measurements no violation of the fluctuation dissipation rela- 
tion is observed. The reasons of this large difference between the electrical and 
rheological measurements are discussed. 



: 1 Introduction 

• Many physical systems in nature are not in thermodynamic equilibrium because they 
. present very slow relaxation processes. A typical example of this phenomenon is the 
aging of glassy materials: when they are quenched from above their glass transition 
temperature Tg to a temperature T < Tg, any response function of these systems de- 
pends on the aging time ta spent at T. For example, the dielectric and elastic constants 
of polymers continue to evolve several years after the quench . Because of these slow 
relaxation processes, the glass is out of equilibrium, and usual thermodynamics does 
not apply. However, as this time evolution is slow, some concepts of the classical ap- 
proach may be useful for understanding the glass aging properties. A widely studied 
question, is how the temperature of these systems can be defined. One possible answer 
comes from the study of the deviation to the Fluctuation Dissipation Relation(FDR) 



in an out of equilibrium system (for a review see ref. p|, ^, ^). In this letter we show 
that this approach is relevant for the study of a colloidal glass formation, where a 
strong violation of FDR is measured. Implications of this observation go beyond the 
physics interest. Indeed FDR is used as a tool to extract, from fluctuations measure- 
ments, several properties in biological, chemical and physical systems 0]. Our 
results indicate that before extending this smart technique to other systems, one has 
to carefully ensure that these systems are in equilibrium. 

In order to understand this new definition of temperature, we have to recall the 
main consequences of FDR in a system which is in thermodynamic equilibrium. We 
consider an observable V of such a system and its conjugate variables q . The response 
function xvqi^), at frequency u = 00/211, describes the variation 6V{u) of V induced by 
a perturbation 6q{Lj) of q, that is Xvqi^) = SV{ij)/6q{uj). FDR relates the fluctuation 
spectral density of V to the response function xvq and the temperature T of the system: 



SiLu) = ^^Im[xv,{cu)] (1) 

TTUJ 

where ^(ti;) =< [^(co')^ > is the fluctuation spectral density of V, ks is the Boltz- 
mann constant, Imlxvqi^)] is the imaginary part of xvqi^)- Textbook examples of 
FDR are Nyquist's formula relating the voltage noise to the electrical resistance and 
the Einstein's relation for Brownian motion relating the particle diffusion coefficient to 
the fluid viscosity 0. 

When the system is not in equilibrium FDR, that is eq.|l|, may fail. For example 
violations, of about a factor of 2, of eq.|I| have been observed in the density fluctuations 
of polymers in the glassy phase 0. The first to propose that the study of the FDR 
violations are relevant for glassy systems was Sompolinsky |]TU|. This idea, which was 



generalized in the context of weak turbulence has been recently reconsidered by 
Cugliandolo and Kurchan and successively tested in many analytical and numerical 
models of glass dynamics 0],[|13|-[|19|. Let us briefly recall the main and general findings 
of these models. Because of the slow dependence on ta of the response functions, it has 
been proposed that eq.|I| can be used to define an effective temperature of the system, 
specifically: 



Im[xvq{ta,uj)\ 

It is clear that if eq.|I| is satisfied Te// = T, otherwise T^ff turns out to be a decreas- 
ing function of and uj. The physical meaning of eq.^ is that there is a time scale 
(for example ta) which allows to separate the fast processes from the slow ones. In 
other words the low frequency modes relax towards the equilibrium value much slower 
than the high frequency ones which rapidly relax to the temperature of the thermal 
bath. Therefore it is conceivable that the slow frequency modes keep memory of higher 
temperatures for a long time and for this reason their temperature should be higher 
than that of the high frequency ones. This striking behavior has been observed in sev- 
eral numerical models of aging 0],[|13|-||19[. Further analytical and numerical studies 
of simple models show that eq.^ is a good definition of temperature in the thermody- 
namic sense P, H- In spite of the large amount of theoretical studies there are only 



a few experiments which show a violation of FDR in real materials p|, 0. However 
these measurements are done at a single frequency and there is no idea on how the 
temperature relaxes as function of time and frequency. The experimental analysis of 
the dependence of Teff{uj,ta) on u and ta is very useful to distinguish among different 



models of aging: FDR violations are model dependent P,||T3|-|T9| . 

For these reasons we have experimentally studied the violation of eq.|l] during the 
colloidal glass formation in Laponite RD a synthetic clay consisting of discoid 
charged particles. It disperses rapidly in water and solidifies even for very low mass 
fraction. Physical properties of this preparation evolves for a long time, even after 
the sol-gel transition, and have shown many similarities with standard glass aging 
22| , p3| . Recent experiments [^] have proved that the structure function of Laponite 



at low concentration (less than 3% mass fraction) is close to that of a glas. As in our 
experiment the Laponite concentration is low, we call the solid like Laponite solution 
either a colloidal glass or simply a glass. 

In our experiment we measure the time evolution of the Laponite electrical and 
rheological properties during the colloidal glass formation. The paper is organized as 
follows. In the next section the electrical measurements are described. Rheological 
measurements are discussed in section 3. Finally we conclude in section 4. 



2 Electrical measurements 
2.1 The experimental apparatus 

The Laponite solution is used as a conductive liquid between the two golden coated 
electrodes of a cell (see fig.l). The Laponite solution is prepared in a clean N2 at- 
mosphere to avoid CO2 and O2 contamination, which perturbs the electrical measure- 
ments. Laponite particles are dissolved at a concentration of 2.5% mass fraction in 
pure water under vigorous stirring during 20mm. To avoid the existence of any initial 
structure in the sol, we pass the solution through a 1/im filter when filling our cell. 
This instant defines the origin of the aging time ta (the filling of the cell takes roughly 
two minutes, which can be considered the maximum inaccuracy of ta)- The sample 
is then sealed so that no pollution or evaporation of the solvent can occur. At this 
concentration, the light scattering experiments show that Laponite structure functions 
are still evolving 500/i after the preparation [^. We only study the beginning of this 
glass formation process. 

The two electrodes of the cell are connected to our measurement system, where we 
alternately record the cell electrical impedance Z(ta, uj) and the voltage noise density 
Sz{ta,uj) (see Fig.|T]). Taking into account that in this configuration Im [xvq{ta,^)] = 
uRe [Z(ta,u!)], one obtains from eq.^ that the effective temperature of the Laponite 
solution as a function of the aging time and frequency is: 



Teff{ta,Uj) = nSzita,Uj)/2kBRe [Z{ta,Uj)] 

which is an extension of the Nyquist formula. 
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Figure 1: Experimental set-up The impedance under test Z is a cell (diameter 
7cm, thickness 3cm) filled with a 2.5wt% Laponite sol. The electrodes of the cell are 
golden coated to avoid oxidation. One of the two electrodes is grounded whereas 
the other is connected to the entrance of a low noise voltage amplifier character- 
ized by a voltage amplification G. With a spectrum analyzer, we alternately record 
the frequency response FR{uj) =< yo^^ jy^^^'^ > (switch A closed) and the spectrum 
S{ijj) =< > (switch A opened). The input voltage V^^'^ is a white noise ex- 

citation, thus from FR{uj) we derive the impedance Z{u)) as a function of u, that is 
Z{u!) = Rq/{G/FR{(jj) — 1); whereas from S{ijj), we can estimate the voltage noise of 
Z, specifically Szi^j) = [S{uj) — Sa{uj)]/G'^ where Sa{u!) is the noise spectral density of 
the amplifier 



2.2 Experimental results 

In Fig.^a), we plot the real and imaginary part of the impedance as a function of 
the frequency u, for a typical experiment. The response of the sample is the sum of 2 
effects: the bulk is purely conductive, the ions of the solution follow the forcing field, 
whereas the interfaces between the solution and the electrodes give mainly a capacitive 
effect due to the presence of the Debye layer This behaviour has been validated 
using a four-electrode potentiostatic technique |^ to make sure that the capacitive 
effect is only due to the surface. In order to test only bulk properties, the geometry 
of the cell is tuned to push the surface contribution to low frequencies: the cutoff 
frequency of the equivalent R-C circuit is less than 0.6Hz. The time evolution of the 
resistance of one of our sample is plotted in Fig.^(b): it is still decaying in a non trivial 
way after 24h, showing that the sample has not reached any equilibrium yet. This 
aging is consistent with that observed in light scattering experiments . 

As the dissipative part of the impedance Re{Z) is weakly time and frequency de- 
pendent, one would expect from the Nyquist formula that so does the voltage noise 
density Sz- But as shown in Fig.||, FDR must be strongly violated for the lowest 
frequencies and earliest times of our experiment: Sz changes by several orders of mag- 
nitude between highest values and the high frequency tail ||26|. This violation is clearly 
illustrated by the behavior of the effective temperature in Fig.H. For long times and 




Figure 2: The response function (a) Frequency dependence of a sample impedance 
for 2 different aging time: ta = 0.3h, real (A) and imaginary (a) part; ta = 24/i, real (•) 
and imaginary (o) part. The equivalent circuit for the cell impedance is shown in the 
inset: Z is the sum of a resistive volume Rb and a capacitive interface Cg between the 
Laponite solution and the electrodes. The increase of Re{Z) toward small frequencies 
< Vc) is due to the dissipative part of the capacitance (loss tangent ~ 0.2). For 
V > Vf. the impedance of the cell is dominated by the bulk resistance Rb- (b) Time 
evolution of the bulk resistance. This long time evolution is the signature of the aging 
of the sol. In spite of the decreasing mobility of Laponite particles in solution during 
the gelation, the electrical conductivity increases. 



high frequencies, the FDR holds and the measured temperature is the room one (SOO-ft'); 
whereas for early times Tgjj climbs up to 3.10^7^' at IHz. Notice that the violation 
extends to frequencies much larger than the {).QHz cutoff separating the volume from 
the surface effects. Moreover, the scaling presented in inset of Fig.|^ seems to indicate 
that Tf.ff can be even larger for lower frequencies and lower aging times. Indeed, we 
found in all the tested samples no evidence of a saturation of this effective temperature 
in our measurement range. In order to be sure that the observed violation is not due to 
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Figure 3: Fluctuations Voltage noise density of one sample for different aging times. 
The strong increase of Sz for low frequencies is quite well fitted by a power law cu", 
with a = —3.3 ± 0.4. This effect is a decreasing function of time, and a good rescaling 
of the data with a cutf law can be achieved as shown in the inset for /5 = 0.5 ± 0.1. 

an artifact of the experimental procedure, we filled the cell with an electrolyte solution 
with pB. close to that of the Laponite sol such that the electrical impedance of the 
cell was the same. Specifically we filled the cell with NaOH solution in water at a 
concentration of 10~^ mol.l~^ . The results of the measurements of Tejj are shown in 
fig.|^ at two different time after the sample preparation. In this case we did not observe 
any violation of FDR at any time. 




Figure 4: Effective temperature of Laponite Effective temperature as a function 
of frequency for different aging times We restrict tlie frequency range to IHz 

to limit tlie surface contribution to tlie measurement. As Sz in Fig.^, T^jf strongly 
increases and reaches huge values for low frequencies and short aging times. This large 
violation is observed in numerical simulations of systems presenting domain growth 
process. 




Figure 5: Effective temperature of an NaOH solution in water. The effective 
temperature is plotted as a function of frequency for two different times after the 
preparation. This solution has a pB. close to that of the Laponite, but no violation is 
observed in this case for any aging time. 



2.3 Discussion 



Let us now briefly discussed the results. The observed very large value of T^jf is 
of course very striking. However the existence of infinite T^ff was predicted ^ and 
numerically verified [jl6[ in systems presenting domain growth process. This is probably 
the way in which the Laponite solution makes the transition towards the colloidal glass 
state. The influence of these domains on the electrical conductivity of Laponite is 
related to the electrostatic interaction among the Laponite disks, which is a widely 
studied and not yet completely understood problem[^. The physical origin of the large 
Tgff certainly lays in these complex interactions. However one may wonder whether the 
observed T^ff is mainly due to bulk effects or to a conductive phenomenon produced 
at the Debye layers on the electrodes. In our opinion the latter has to be excluded for 
two reasons. The first one is that the violation extends to frequencies much larger than 
the 0.6Hz cutoff separating the volume from the surface effects. The second reason is 
that complementary measurements, done in cells with different geometries, show that 
for frequency larger than the cutoff Uc the effect of the Debye layers on the violation is 
negligible. 

The comparison with the results of the other experiments ||^, ^] is difficult because 
these two experiments are done in different materials and at a single frequency. Thus 
in these two experiments it is impossible to know the evolution of the violation both 
in frequency and time. Furthermore the violation observed in these experiments is of 
order one whereas in our experiment it is several order of magnitudes. However our 
experiment and the one on supercooled fluid ||2^ have a common important result, 
which merits to be stressed. The violation of FDR at frequency u persists for times 
which are several order of magnitude larger then the time l/w, which is often considered 
as the characteristic time of the violation of FDR at frequency uj. 



3 Rheological Measurements 

In the previous section we have shown that the FDR is strongly violated by the electrical 
properties of Laponite. We want to understand whether a violation can be observed 
in the measurements of other physical properties. Of course there are no reasons to 
assume that T^ff is the same for all the variables but the differences and the analogies 
found in the evolution of T^ff obtained by the measurements of several variables may 
give new insight on aging theories and on the coupling of the different variables in an 
aging material. Therefore we performed rheological measurements on Laponite and we 
checked FDR in these measurements. 



3.1 Experimental apparatus 

To achieve this result we built a new rheometer which is sensitive to thermal fluctua- 
tions. The principle of the rheological measurement is a standard one and is illustrated 
in fig.p. We describe here only the main features, more details can be found in ref.[29|. 
A rotor of diameter 12mm is inserted in a cylindrical cell. The gap of 1mm between the 
rotor surface and the cell is filled with the fluid under study. The rotor is suspended by 
two steel wires. On the top of the rotor we fixed an optical prism. This prism is part 



of a Nomarski interferometer which is used to measure the rotation angle 9 of the 
rotor. The sensitivity of this system is better than 10~^^rad/\/Hz corresponding to a 
torque on the rotor of about 10~^^N ■ m/ V Hz. An external torque Text can be applied 
to the rotor by using the electrostatic interaction of a capacitor (see fig.P). One of the 
two electrodes of the capacitor is fixed on the rotor whereas the other is fixed on the 
cell walls. When a voltage difference Vc is applied on this capacitor the attractive force 
between these two electrodes produces a torque on the rotor which is equilibrated by 
the stiffness of the steel wires. 




filling hole 



a) 



b) 



Figure 6: Rheometer Vertical (a) and horizontal (b) cross section of the rheometer. 
The torque is applied through the electrostatic interaction of electrodes Ai_2 and Bi_2 
when the voltage Vc is applied to this capacitor. 



Let us consider the simple case of a Newtonian fluid of viscosity rj. The response 
of this torsion pendulum in Fourier space is 



X,dText 
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(4) 



where J is the rotor inertia moment, k the steel wires stiffness and a is a geometric 
factor. 

In the experiment we have access to two quantities : the angular position 9 of the 
rotor, and the voltage Vc driving the torque. Text being quadratic in Vc, we add a 
constant offset to it in order to linearise this relation. For small displacements, the 
response function XsText is then simply proportional to XeVc- The missing constant can 
be found by performing a inertial calibration of the response : the real part of 1/ XeVext 
is a parabola whose quadratic coefficient is the rotor inertia moment J. J is known 
with a good precision, and thus can be used to calibrate the measurement. 



This rheometer has been tested using a sihcon oil with rj = 2Pa s. We first measure 
the response by using a white noise voltage exitation for Vc with an amplitude corre- 
sponding to a torque Tf^^t — 10~^^N m/ \/Hz. The real and imaginary part of '^/xeVext 
are plotted as a function of frequency in fig.^(a) and fig.|^(b) respectively. We see that 
in agreement with eq.^ the real part of 1/ XeVe^t is very well fitted by a parabola whereas 
the imaginary part is linear for small v. The deviation of the data from the linear be- 
haviour is due to the non-newtonian character of the silicon oil at higher frequencies. 
Once the response is known one can set the external torque to zero and measure Sq, 
the spectrum of the thermal fluctuations of 6. In this case the fluctuation dissipation 
relation is : 

Se = —^Irn{xev.J (5) 

UJ 

By inserting the measured Xev^^t 6q.|^, we get an estimation of the fluctuation 
spectrum Sq. The comparison between this computed Sq and the measured one is 
done in fig.H. Except for the existence of the peaks due to environmental noise the 
agreement is quite good. Thus the rheometer has enough sensitivity to verify FDR in 
viscous fluids. 




Figure 7: Response of the rheometer to a white noise excitation. The real (a) and 
imaginary (b) part of '^-/xeVext plotted as a function of frequency. The rheometer 
is filled with a viscous oil with rj = 2Pa s. The dashed curve in (a) corresponds to a 
quadratic fit of the data, and allows inertial calibration of the measurement (see text 
for details). The dashed line in (b) is a linear fit of low frequency data, it is consistent 
with a newtonian behavior of the fluid. 



3.2 Measure on Laponite 

We have studied the rheological properties of Laponite at 3% mass concentration in 
water. Laponite has been prepared in the way described in section 2) for electrical 
measurements and the rheometer is mounted inside a container filled with a clean 
Nitrogen atmosphere. Moreover, a thin oil layer on top of the Laponite solution ensures 
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Figure 8: Spectrum of the thermal fluctuations of Q. The rheometer is filled with the 
same oil used to measure the response reported in fig.^.The continuous line is the result 
of the measurement whereas the dashed line is computed inserting in eq.^ the measured 
response of the rheometer and T = SOOii'. 

no evaporation can occur. We first measure the response at different times after the 
preparation. The measured elastic and viscous modules of Laponite are plotted as a 
function of frequency in fig.^ (a) and (b) respectively. We clearly see that at short 
the liquid like solution has only viscous response. As time goes on an elastic modulus 
appears and the viscosity grows of about one order of magnitude. After 24 hours the 
solution has a solid like aspect with a viscoelastic response. The measurement the 
thermal fluctuation spectrum 5*^ averaged on the first hour of the aging processes of 
Laponite is plotted in fig.^. The peak at X.'bHz has no physical meaning and is due 
to a mechanical resonance of the table. Whereas the resonance at IHz is the resonant 
frequency of the torsion pendulum of the rheometer. The measure is compared to 
the prediction of the fluctuation dissipation theorem by inserting in eq.^ the measured 
values of the Laponite viscous response and T = 300i^. We also inserted in eq.| the 
effective temperature obtained by the electrical measurements averaged on the first 
hour. We clearly see that in the case of rheological measurements no violation of FDR 
can be detected within experimental errors. If a violation exists, it is much smaller 
than that observed in the electric measurements. 

4 Discussion and Conclusion 

We have studied the fluctuation dissipation relations during the transition of Laponite 
from a fluid-like state to a solid-like colloidal glass using electric and rheological mea- 
surements. Strong aging properties have been observed both in the electrical and 
rheological response functions. The behaviour of thermal fluctuations is instead quite 
different. In electric measurements a large violation of the dissipation fluctuation re- 
lation is observed at the beginning of the transition from a fluid like sol to a glass. As 
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Figure 9: Aging of the viscoelastic module of Laponite at a 3% mass concentration in 
water. The real (a) and imaginary (b) part of the viscoelastic module are plotted as a 
function of frequency. 
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Figure 10: Thermal fluctuation spectrum of the rheometer filled with Laponite. The 
line with o is the the result of the measurement. The line with x is the FDT prediction 
for T = 300K. The line with + is obtained from eq.^ by inserting in T the values of 
Teff estimated from electric measurements. 



predicted by the theory the amplitude and the persistence time of this violation are 
decreasing functions of frequency. The observed effective temperature reaches lO^i^' 
for frequencies smaller than IHz. This violation involves frequencies which are much 
larger than l/tw Instead in the rheological measurements no detectable violation of 



the fluctuation dissipation relation has been observed. The reasons of this behaviour 
and of such a big difference between the two physical properties are not understood. 
The first striking result is the very high effective temperature measured in the electrical 
properties. These very high effective temperatures has been observed in the numerical 
simulation of aging models presenting domain growth. Thus this behaviour of electrical 
properties of laponite could be explained by the formation of large domains during the 
transition toward the colloidal glass state. However one wonders whether the growth 
of these domains does not affect the fluctuations in the rheological measurements. An- 
other possible explanation of the large Tg// in electric measurements is the dissolution 
of aggregated particles and consequently of ions in the Laponite-water solution. It 
is conceivable that this effect should not influence the rheological measurements. The 
second striking result, which merits to be discussed, is the fact that in electric measure- 
ments the violation involves frequencies which are much larger than 1/tw This results 
agrees with the experiment of ref. but it is different from what is theoretically pre- 
dicted and computed in numerical simulation. This observation leads us to another 
possible explanation of the discrepancy between electrical and rheological properties 
of Laponite. Indeed the time scale of relaxation of the fluctuations could not be the 
same for the two properties. At the moment we have not a clear and precise answer 
to this problem and more precise measurements are necessary in order to understand 
this large difference between rheological and electrical measurements. 

The results of this paper are certainly preliminary but they show that the measure- 
ments of FDR in aging systems can give more insight to the problem of material aging. 
Indeed aging has been always characterized only by the measurements of the response 
function. We have shown here that the association of this measurement to that of the 
thermal fluctuations may show new and unattended aspects of an aging material. 
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